CUBATURE FORMULAE FOR ORTHOGONAL POLYNOMIALS 

IN TERMS OF 
ELEMENTS OF FINITE ORDER OF COMPACT SIMPLE LIE GROUPS 

ROBERT V. MOODY AND JIRI PATERA 

Abstract. The paper contains a generalization of known properties of Chebyshev polyno- 
mials of the second kind in one variable to polynomials of n variables based on the root 
lattices of compact simple Lie groups G of any type and of rank n. The results, inspired by 
work of H. Li and Y. Xu where they derived cubature formulae from ^4-type lattices, yield 
Gaussian cubature formulae for each simple Lie group G based on interpolation points that 
arise from regular elements of finite order in G. The polynomials arise from the irreducible 
characters of G and the interpolation points as common zeros of certain finite subsets of 
these characters. The consistent use of Lie theoretical methods reveals the central ideas 
clearly and allows for a simple uniform development of the subject. Furthermore it points 
to genuine and perhaps far reaching Lie theoretical connections. 
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1. Introduction 



During most of the century and half long history of Chebyshev polynomials, only polyno- 
mials of one variable were studied |17j . In recent years a considerable overlap of the subject 
,__! can be found in the emerging and practically important field of cosine and sine transforms 

J> [16. 18J. In the absence of additional constraints, truly higher dimensional generalizations of 

Chebyshev polynomials are hidden in the vast number of possibilities of defining orthogonal 
■^ polynomials of more than one variable |19l [10] . In 2D [9] such constraints were provided by 

fNJ requiring that the polynomials of two variables be simultaneous eigenvectors of two differential 

iy-^ operators. 

^^ Our motivation here comes from two directions: 

(i) The method of constructing Chebyshev-like polynomials for any simple Lie group [H] , 
and particularly from the recognition of the basic role played by the group characters. 
.£h (ii) The work of H. Li and Y. Xu in which they derived cubature formulae based on the 

/\ symmetries of ^4-type lattices [8J . 

From p3] we understood that a general formulation, uniform over all the types and ranks of 
simple Lie algebras should be possible, and from [8] we saw what possibilities, beyond the 
construction of the polynomials, might be achievable in a general formulation guided by the 
theory of compact simple Lie groups. 

In this paper the characters of the irreducible representations of a compact simply-connected 
simple Lie group G of rank n play the role of the Chebyshev polynomials (of the second kind) 
and elements of finite order |llj in G give rise to interpolation points through which we arrive 
at Gaussian cubature formulas. The ring generated by the characters of G has a Z-basis con- 
sisting of the irreducible characters and it is a polynomial ring in terms of the n irreducible 
characters of the fundamental representations. These fundamental characters then serve as 
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new variables for functions denned on a bounded domain fi C M n . This domain fi is derived 
from the fundamental domain of the affine Weyl group of G and the kernel function used is 
the absolute value of the denominator of Weyl's character formula [20] . 

There are two technical ingredients in the paper, which are indispensable for the uniformity 
of our approach to simple Lie groups of all types and thus for generality of our conclusions. 

(i) The 'natural' grading of polynomials by their total degree is replaced by new m-degree 
grading. It is based on a set of Lie theoretical invariants, called the marks, which are unique 
for simple Lie group of each type. The two gradings coincide only in the case of A n . 

(ii) The set of interpolation points is uniformly specified for each simple Lie group as a finite 
set of lattice points, characterising all conjugacy classes of elements of certain order in the 
underlying Lie group. 



The cubature formula (see Cor. 7.2) is then a formula that equates an integral of a general 
polynomial P, say of m-degree < 2M+1, to weighted finite sums of the values of P sampled at 
the interpolation points which are common zeros of the polynomials of m-degree M + 1. The 
cubature is Gaussian [21] in the sense that the number of the interpolation points coincides 
with the dimension of the space of polynomials of m-degree < M. It turns out that the 
interpolation points are the elements whose adjoint order divides M + h where h is the 
Coxeter number. 

In more detail, having fixed a specific m-degree M, we are interested in the properties of 
the set of all polynomials of m-degree at most M. There are three main results: 



1) Thm. 6.1 The interpolation points are common zeros for the set of polynomials of 
m-degree M + 1. These interpolation points are directly related to regular elements 
of finite adjoint order M + h in G, where h is the Coxeter number of G, and their 
number is precisely the dimension of the space of polynomials of m-degree less than 
or e qual to M. 



2) Cor. 7.2 There is a cubature formula that equates weighted integrals J^fK 1 ' 2 of 
each polynomial / of m-degree < 2M + 1 with if -weighted linear combinations of its 
values at the interpolation points. 



3) Prop. 8.1 The expansions of functions via irreducible characters (or rather these 
characters interpreted as polynomials) using interpolation points yield the best ap- 
proximations in the K 1 > 2 -weighted L 2 -norm on the functions on f2 . 

These results and their proofs are natural and completely uniform within the framework 
of the character theory of simple Lie groups and their elements of finite order. In fact it is 
this naturalness and perfection of fit that suggests that there are a deeper Lie theoretical 
implications to all of this that are still to be discovered. The potential role of the finite 
reflection groups in the theory of orthogonal polynomials has been recognized already in 
|3j, although only a limited use of the groups is made there, the objects of interest being 
group invariant differential and difference-differential operators. The present paper can be 
understood as a new contribution to the fulfilment of that potential, much closer to the 
properties of the simple Lie groups which give rise to the finite reflection groups. 

Our discussion requires a certain familiarity with root and weight lattices, their correspond- 
ing co-root and co-weight lattices, the Weyl group and its affine extension, and particularly 
the structure of the natural fundamental domain for the affine Weyl group. We use the £|2j 
which sets up the notation, to briefly review the key points of this theory, while £J3] and $4] 
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contain preparatory extensions of the standard theory. Consideration of the subject of the 
paper really starts in £J5j 

2. Simple compact Lie groups 

This section is a brief review of the material that we need for this paper and establishes 
the notation that we shall be using. The main facts about simple Lie groups and their 
representations are classical and can be found in many places. One source that uses the 
same notation as in this paper is [5], Vol.1. Material on the fundamental domain, and in 
particular information about the stabilizers of its various points, can be found in Ch.V of [1\ 
and material on both the fundamental domain and the elements of finite order can be found 



m 

Let G be a simply connected simple Lie group with Lie algebra q. We let cq denote the 
order of its centre and let T be a maximal torus of G. We let it be the Lie algebra of T, so 
that we have the exact sequence 

o _► _► t "^"H T — ► l 

via the exponential map. Using t instead of the actual Lie algebra it of T has the advantages 
that the Killing form (• | •) of G restricted to t is positive definite and e 2nl ^'' is perfect for 
the Fourier analysis to follow. Let n := dim^t, the rank of G. 

The kernel of e 2m ^' on t is the co-root lattice Q, and t/Q ~ T naturally expresses T as 
a real space factored by a lattice. We denote by t* the dual space of i and let (• , •) be the 
natural pairing of t* and i. 

Let Of be the Haar measure on T that gives it volume equal to 1 . In practice we most often 
write integration over T in the form of integration over some fundamental region FR for Q 
in t: 

(1) I fd9 T = I f{e 2 ™)dOi 

JT JFR 

where 9{ is ordinary Lebesgue measure in t normalized so that FR has volume equal to 1. In 
the sequel no fundamental domain FR ever makes an appearance, but rather we work with 
a smaller fundamental domain F° of the affine Weyl group, see below. 

Given any finite dimensional complex representation V of the Lie group G, the action of 
the elements of T on V can be simultaneously diagonalized, the resulting eigenspaces in V 
being the weight spaces. The corresponding action of t on these weight spaces then affords 
elements A € i* for which x S t acts as exp(27ri(A, x)). Naturally (X,Q) C 7L and the set of 
weights taken over all finite dimensional representations is the subgroup P C t* which is the 
Z-dual of the co-root lattice relative to our pairing. P is the weight lattice of G relative to T. 
The adjoint representation, the representation of G on its own Lie algebra, produces its own 
set of weights. Apart from the weight 0, which is of multiplicity n, the remaining weights 
of the adjoint representation are the roots of G and they all occur with multiplicity 1. They 
generate the root lattice Q of G, and we have Q C P with index equal to cq defined above. 
The Z-dual of Q in t is the co-weight lattice P D Q with the index of Q in P being cq, again. 

Let W be the Weyl group associated with T, namely the quotient of the normalizer of T 
in G. Then W acts as a finite group of isometries of i relative to (• | •), and W stabilizes Q. 
It also stabilizes P, and if we pull its action over to t* by duality, then W also stabilizes the 
weight lattice P and the root lattice Q. 
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The relationships between the lattices and between the various root and weight bases and 
their co-equivalents described below are summarized in: 

{«!,..., a n } C 
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n x n 
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t* 


t 



(2) P P D {ut,...,0J n } 



The times symbol is meant to indicate that Q and P, as well as Q and P, are in Z-duality 
with each other. The indicated bases are also in Z-duality. 

We have the semi-direct product W & q = W ix Q acting on t, with Q acting as translations 
and the Weyl group as point symmetries. W & q is called the affine Weyl group. 

A fundamental region F for t under the action of W & g can be given as follows. Let A be 
the set of all the roots and let IT = {ai, . . . , a n } be a simple system of roots for A. Let A+ 
(resp. A_) denote the corresponding set of positive (resp. negative) roots, and let 

(3) a = -(m\a\ H \- m n a n ) 

be the lowest rootjj of A. The positive integers mi , . . . , m n are called the marks of G. The 
marks and co-marks are shown on Fig. [TJ They are independent of the choice of T and the 
choice of II (but depend on the choice of ordering of a\, ■ ■ ■ ,a n ). It is convenient to define 
one more mark, ttlq := 1, so one has YJj=o m j a j = 0- One then has the Coxeter number of G 

n 

j=0 

This constant appears in many places and in many guises in the theory, and plays an important 
role in what follows. 
Now we define 

(4) F := {x e 1 1 (aj,x) > for all j = l,...,n, {-a ,x} < 1} . 

This is a simplex in t. The Weyl group W contains the so-called simple reflections n , . . . , r n 
in the walls 

(5) Hj :={xet\ (aj,x) = 0; j = l,...,n} 

and W is generated by these reflections. The length function on W is the homomorphism 

l:W — > {±1}, rj ^ -1 for all j = 1, . . . ,n . 

There is a unique element u>opp £ W which maps A + into A_. It is an involution and is 
minimally represented by the product of exactly | A + | of the simple reflections. 

W a Q is also a reflection group and is obtained by adding to n, . . . ,r n the additional gen- 
erator ro which is reflection in the wall 

(6) H :={xei\(-a ,x) = l}. 

There is a similar length function for W & q The action of W & ^ on F tiles the entire space t 
with copies of F, and in this way F serves as a fundamental domain for it. See Fig. [2] for an 
example. 



In what follows, the use of the lowest root rather than the highest root, which might seem more natural, 
is in keeping with notation from the theory of affine root systems. 



CUBATURE FORMULAE FOR SIMPLE LIE GROUPS 5 

The way in which F is the fundamental region is rather beautiful: 

• Each Wog orbit in t has a unique element in F. 

• For x £ F the stabilizer of x in W a g is the subgroup of W & g generated by the reflections 
r ji i = 0, 1, . . . , n, for which x is on the wall Hj. In particular, all x £ F° have trivial 
stabilizer in W a g. 

The fundamental co-weights are the elements w^Gt dual to the simple roots: (aij, W&) = <5jfc. 
They form a Z-basis of P. In terms of them we can write F as the convex hull of 

mi rri2 m ri 

Although each element of W and each element of W & g is a product of reflections, it is not 
uniquely so. However, the parity of the number of reflections is well-defined and we use the 
symbol (— l)" w > to indicate this parity, it being ±1 according to evenness or oddness. 

The Cartan matrix of G is the integer n x n matrix 

. . . ( %OLi I (Xj)\ 

A=(A ij )= v ' J M , l<i,j<n. 

V ("J I OLj) J 

It is unique up to the numbering of simple roots. The Cartan matrices classify the compact 
simply-connected simple Lie groups into the well-known A, B, C, D series and the five 
exceptional groups Eq, Ej, Eg, P4, and G2. 

The simple co-roots {di, . . . , a n } C Q corresponding to {a±, . . . , a n } C Q are defined by 

2(«j \ dj) 
{ai,aj} = A^ = : — J — tor all 1,3 . 

{OLj I Oj) 

They form a Z-basis of Q and their ^-translates in Q form the set of co-roots A. Actually 
A and the simple co-roots are a root system and simple roots for a simply-connected group 
G whose Cartan matrix is A . We do not need this group directly in what follows but we 
occasionally use information about co-objects that we know is true from the fact that they 
have such an interpretation. 

Dual to the co-roots we have the fundamental weights coj £ P defined by (cjj,ak) = Sjk- 
The fundamental weights form a Z-basis of P. 

Each finite dimensional irreducible representation L of G has a unique one dimensional 
weight space L x , A £ P, with the property that all other weights of L are of the form A — /3, 
where /3 is a sum of positive roots. We have 

A G P+ := {fi g t* I {n, on) S Z^°, i = 1, . . . , n} . 

Here P + C P is the set of dominant weights. The dominant weight A is called the highest 
weight of L. If we designate L now by L(\), then the correspondence 

AeP + ^ L(\) 

classifies all the irreducible finite dimensional representations of G up to isomorphism. 
We shall also need the set of strictly dominant weights 

P++ :={(iet* I ( fi ,a i )€Z >0 , t = l,...,n}. 
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Figure 1. The usual Coxeter-Dynkin diagrams, which are simply a visually 
efficient way to encode the Cartan matrices, are shown, along with information 
giving the corresponding marks and co-marks. The circular nodes, ignoring 
those with a dot in them, stand for the simple roots with the convention that 
open (resp. filled) circles indicate long (resp. short) roots. When there are 
both types of nodes, interchanging open and filled results in the co-diagram. 
Thus types B and C interchange and types F and G end up as the same type 
but with the numbering of the roots permuted. The dotted node stands for the 
root ao (which is linearly dependent on the simple roots with the negatives of 
the marks as its coefficients). Under the duality operation of roots to co-roots, 
it passes from being the lowest long root to being the lowest short co-root, a 
fact that we do not explicitly have to use here. 

The links between roots occur only when the roots are not orthogonal to one 
another. The marks and co-marks are shown as the fractions m/rh attached 
to the corresponding nodes of the diagram. When both are equal to one, they 
are not shown. We refer the reader to [2] for more details. The numbering of 
simple roots goes from the left to right, the node above the main line carrying 
the highest value. Dotted cicle has number 0. 



The 'simplest' strictly dominant weight is p defined by (p, &j) = 1 for i = l,...,n. The 
element p plays an important role in what follows. We also know that 
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The character of a finite dimensional representation L of G is the mapping 

9 •-> tr L g , g e G. 
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Since the character is unaffected by conjugation by group elements, we can always restrict it 
to T without any loss of information. We then further consider it as a function XL on t by 

XL- x ^ tr L exp(27rix) = ^dimLV^' x) , 

where /U runs over the weights of L. In particular, we have the characters x\ = Xl(\) > ^ £ P + ■ 
Weyl's character formula is 

() XA( )_ £^(-i)W™<^> "• W ' fora11 xeu 

Also there is the special product formula for the denominator: 

(9) S p {x) = Y^ (-lf w h 2lTi{wp ' x) = Yl U ni{a ' x) - e ~ m{a ^ 

The functions S\ +p (x) and S p (x) are W^-skew invariant whereas their quotient is W- 
invariant. The functions S\ +P (x) = ^2 weW (—l) l ^e 2m ^ x+p ^' x ^ are called 5-functions in 
[T3"l [H] due to their similarity in form to the sine function, which they are in the case 
of A\. In any case the values of the characters and of the ^-functions are determined by 
their values on the fundamental domain F, and it is this fact that becomes the centre of our 
attention in the sequel. 

Note that none of the reflecting hyperplanes ([5j or ([6]), nor indeed any reflecting hyperplane 
of W a g, meets the interior F° of F and so S p (x) is never in F°. Thus IS^x)) 2 is positive 
on the interior F° of F. On the other hand, S\ +p (x) and S p (x) vanish on the boundary of F 
since x £ Hj implies that rjX = x while replacing x by VjX in any S-function changes its sign. 

3. VF-INVARIANT AND VF-SKEW INVARIANT FUNCTIONS ON T 

3.1. The algebra of formal exponentials. 

Starting with the weight lattice P one may form the algebra C[P] of formal exponentials, 
which has a C-basis of symbols e A , A G P, together with a multiplication defined by bilinear 
extension of the rule 

e A e M = e A+^ _ 

Thus typical elements of C[P] are finite complex linear combinations ^ c\e . C[P] is an 
unique factorization domain and its group of invertible elements are the elements ce A , where 
c /0,Ag P. 

W acts on P and hence as a linear operator and even an automorphism on C[P] by w.e = 
e . There is a partial order on P with A > \i if and only if /x = A — /3 where /3 is a 
sum (possibly empty ) of positive roots. This is important because in the weight systems of 
irreducible representations of G the highest weight is highest in this sense. We will use this 
notion of highest below. 

There are some advantages to introducing formal exponentials at times since they often 
clarify the mathematics. However, we are really interested in their manifestations as functions 
on t and on T which arise from 

(10) x h+ e 2nix ^ e 2ni{x ' x) . 

These mappings are the characters <fi\ : x t— > exp(2iriX.x) of the torus, so C[P], as functions, 
is its character ring. This relates directly back to the previous section where we have defined 
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the G-characters and S- functions, all of which may be viewed as arising from corresponding 
elements Xx, S\ +p of C[P] . 

At the level of functions the action of W is completely consistent: 

w.<j)\(x) = 4>\(w~ x) = exp(27rz X.w~ x) = exp27rzuiA.x = 4> w \(x) . 

We are most interested in the subring C[P] of TV-invariant elements. The simplest forms 
of TV- invariant functions are the orbit sums ^lugw 4> w \, called C- functions in [151 El 031 HI] 
due to their similarity in form to the cosine function (which they are in type A\). More 
relevant here are the characters \\ °f G (restricted to T) already introduced in q2j 

(ii) X \ = Y, dimL ( x ^ e2wi>M 

fieP 

where L(X) is the irreducible representation of highest weight A E P + and dimL(A) M is the 
dimension of its //-weight space. These are VF-invariant and they form a basis for C[P] W . 

Proposition 3.1. [lJ(Ch.VI) CfP] 1 ^ is a polynomial ring with the fundamental characters 
Xwi » • • • Xuj n as the generators Q 

This result really underlies the results of this paper. It says that the fundamental characters 
of G can be used as new variables by which the algebra of invariant elements of C [P] becomes 
a polynomial ring in these variables. It is in working out the Fourier and functional analysis 
implied by this statement that the cubature formulae arise. 

In the sequel we prefer to reserve the word character for the characters of G (as opposed 
to the characters of T) since they are of fundamental importance to the paper. 

3.2. Skew invariants elements of C[P]. 

Elements £ S C[P] for which w.£ = (— ly( w )£ for all w G W are called skew-invariants. 
They play a vital role in the paper. The simplest example is S„, and it is the foundation for 
all the skew-invariant elements. 

Proposition 3.2. [1] (Ch.VI) 5pC[P] w is the set of all W -skew-invariant elements of 

C[P\. 

Later on, when we use the basic characters Xuj- as new variables Xj and have polynomial 
functions of the Xj , we shall have need of a Jacobian for the switch of variables from the Xj 
back to the variables that parameterize T. We establish the key result here. 

For each d£Q there is a unique derivation^ D& on C[P] satisfying 

D & {e x ) = (A,d)e A 

for all A. D^ is a linear in d. 

Let Xuii ; • • • > Xu)„ De the basic characters and let &i, . . . , a n be the standard basis of Q dual 
to {ui, . . . ,uj n }. Let J be the matrix with entries Jjk = D aj Xui k - 

Proposition 3.3. (Steinberaj 

det(J) = S p . 



The orbit sums over the various dominant weights also form a basis for C[P] W and those for u)\, . . . , u> n also 
form a set of generators for it as a polynomial ring. The relationship between orbit sums and corresponding 
characters is a triangular matrix of integers with Is down the diagonal [2]. 

3 A derivation of C[P] is a linear mapping D : C[P] — > C[P] satisfying D {<)>(,) = £>(»£ + <j>D{£) for all 
0,£6C[P]. 

This result is an Exercise to Ch. VI in [l] and attributed to R. Steinberg there. 
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Proof. All the exponentials in Xu) k are of the form e UJk ~^ k where (3^ ls a sum of positive roots, 
and the highest term is e Wj . Now Da.Xu k ^ s a sum OI " exponentials of the same form, but since 
{uj, 6ik) = 8jk> the highest terms only survive along the diagonal of J. Thus when we compute 
the determinant we obtain a sum of signed products e UJl ~^ 1 . . . e^™ - " 1 and only the term from 
the diagonal can contribute an exponential of the form e Wl . . . e Wn , and its coefficient is 1. 
Thus det(J) is an element of C[P] whose highest term is ui\ + ■ ■ ■ + uj n = p and this occurs 
with coefficient equal to 1. 



We shall prove that det(J) is W-skew invariant. Then by Prop. 3.2 it is a multiple of S p . 
Because all the weights in the expansion of det(J) are less than or equal p, this multiple can 
only be a scalar. Since the leading coefficient is 1 in both cases, det( J) = S p , as we wish to 
prove. 

A simple computation shows that wDaW^ 1 = D W a for all a. Fix any I = 1, . . . , n and let 
H = fa r We have natj = a.j — Aij&i. Then, 

ndet(J) = det(r,(J))=det((r J (13a i x Wk ))) 

= det((riD &j ririXui k )) = det((D n a J x^J) 

= dek{(D & .-A lj & l (Xu>k))) = det (( D & 3 ~ A ijD &l ){xu k ))) 

= det((-D dj .(xu,J - AijD &l (Xu h ))) , 

where in the second line we used the W-invariance of the characters. The operation has 
resulted in altering all rows except the Ith row by a multiple of the Zth row (which does not 
alter the determinant) and replacing the Ith row by its negative, which changes the sign of 
the determinant. Thus r\ det(J) = — det(J), which gives the desired skew-symmetry. □ 

3.3. An inner product on C[P] W . 

The natural inner product on C[P] is (■ , -)t defined by 



(f,9h= / fgdOj, 

Jt 

where the 9j is the normalized Haar measure of the torus. Relative to this the functions 4>\ 
form an orthonormal basis. Using this inner product we can complete C[P] in the correspond- 
ing L 2 -norm to the Hilbert space L 2 (T,8j) with the normalized (fi\ forming an orthonormal 
basis, in the sense of Hilbert spaces. Of course we can look at the closure of Cf-P] 1 ^ in 
L 2 (T, 0j), which is in fact the subspace of W^-invariant elements L 2 (T, 9j) of L 2 (T, 9j). 

However, the inner product (• , -)t is not ideal for this subspace, and rather we would like 
to find one with respect to which the characters xa form an orthonormal base. 

We note that for any / G L 2 (T, 6j) , fS„ £ L 2 (T, 6j), and it is skew-invariant with respect 
to W. Form its Fourier expansion 



fS p — 2 ,(fSpi' 



Vt <P» ' 



equality being in the L 2 sense. The summands can be gathered together into VF-orbits, and 
on each orbit the coefficients (fS p ,4> p )j are equal in absolute value and alternate in sign 
according to the parity of the Weyl group elements. The only orbits that do not vanish are 
those contain an weight p G P ++ , and we get 

fSp = 2^ (fS p , 4>x+ P )j S\ +p . 
\eP+ 
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Dividing out the function S„, which is valid as long as the functions are restricted to the 
interior of the fundamental chamber F°, we obtain 

(12) /= ^(/5 p ,0 a+p ) tX a. 

AeP+ 

Now, using the VF-invariance of 6j and the skew-invariance of fS p , we have 



P <Pw(A+p) dOj 



{fS p , 0a +p )t = / fS p 4>x +P Mt = tL I E (-!) W A 

(13) = T^r / / S p Sx+p d8 T = — fxx S P S P d6r= fxx S P S P d9 T . 



Thus (12) and (13) show that / G L 2 (T,9j) w has a Fourier expansion in terms of the 
characters xx, A G P^, with coefficients given by a new inner product defined on L 2 (T,9j) w 
by 

(14) (/,<?) ^ / fgSpSpdBj. 

J F° 



We can then rewrite (12) as 



(15) /= £(/,Xa)xa- 

AeP+ 

We shall use these results in f|8j 

4. Elements of Finite Order 

Elements of finite order [11] (EFOs) in G are used to create the interpolation points for 
the discrete Fourier analysis |12j and cubature formulae to follow. 

We have seen that every element of G is conjugate to one of the form g = exp(27rix) where 
x £ F. The element g is called regular if its centralizer is of dimension n, the rank of G. Since 
each element lies in a torus this is the smallest possible dimension for a centralizer. Regularity 
is a property of the entire conjugacy class of an element and for x £ F it is equivalent to 
saying that x £ F°. 

The condition that g has finite order dividing M is the equivalent to the condition that 
exp(27rzx) = exp(2-7riMx) acts trivially on every irreducible representation, and for this all 
we need is that it acts trivially on every weight space L(X) fl . In turn this requires precisely 
that (/i, Mx) G Z for all weights of P, and finally it is equivalent to Mx G Q, since Q is the 
Z-dual of P. 

In fact what we are going to need here is not that x G -^Q but rather that 

M 

a statement that is equivalent to saying that Ad(g) = 1, i.e. g acts trivially in the adjoint 
representation. In this case we say that g has Ad- order or adjoint order M, even though the 
actual adjoint order, which we shall call the strict adjoint order, namely the least N for which 
Ad(g) = 1 may be some proper divisor N of M. We also say that x G i is an element of 
adjoint order M if exp(27rix) is of adjoint order M. 
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Given the definition above, the conjugacy classes of elements of adjoint order M are rep- 
resented by the points x of the form 

1/ Z^^3 

3=1 
where 

1 n 
(17) Sj G Z- for all j and (— ao,x) = — > rn j s j — 1 

3=1 



1 

(16) x= M^Z s - iLJl 



The regular conjugacy classes of adjoint order M are represented by ( 17 ) where the inequalities 
are made strict. We write Fm (resp. F M ) for the elements of F (resp. F°) of j4d-order M. 
Using rriQ = 1 defined in §|2|we can define so £ Z-° so that 

n 
(18) ^m jSj =M. 

3=0 

Listing all the elements of Fm (resp. F M ) is then just a question of finding all non-negative 



(resp. positive) integer solutions [sq, si, . . . , s n ] to (18). We call [sq, si, . . . , s n ] the ifac coor- 
dinates of x. 

We will be particularly interested in the set Fm+h of elements x € F of j4d-order M + h 
for some non-negative integer M: 

n 

(19) x = M^isx&i + s 2 w 2 H h s n w n ) G jj^-P , where ^J m^ <M + h 

0=1 

where si > 0, . . . , s n > are integers such that X/?=i m j s .?' < Af + ^- Alternatively we have 
the Kac coordinates [so, si, . . . , s n ]. 



Each of the following three conditions assures that x of (19) is in F M+h : 
Sj>0, j = 0,l,...,n; 

n 

2 rn^tj = M , tj := Sj — 1 > 0; 

(20) ,=o 

n 

y m^ij < M (so to completes the sum to M) . 

3=1 

When M = 0, it contains only the element given by so = si = • • • = s n = 1. For M > 0, 
it clearly contains |-F^ +/l | = |i*Af I points. Formulas for the cardinality of \Fm+h\ have been 
worked out for all M and for all simple G in [I]. 

5. Points of -F^+h AS zeros of S-functions 

We are now at a point where we begin the main development of the paper. We fix, once 
and for all a non-negative integer M. The first step is to show that the points of F M+h are 
common zeros of a certain set of S- functions. These points are the interpolation points for 
the cubature formulae to follow. 
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Consider a dominant weight A = A]Wi + • • • + X n co n . We want to find points x £ F^ +h at 
which the S- function S\ +P (x) vanishes: 

g ( x \ — \^ (_l\liw) e 2m{w{\+p),x) _ V^ l_\\l{w) e 2m{\+p,w- 1 x) _ q 

wew wew 

One way to make this happen is to have 

(21) (A + p,x)-{\ + p, rx) £ Z , for all x £ jj^P , 

where r is the reflection in the highest coroot, for if this is the case then the sum collapses in 
pairs adding up to zero. Via (wa,f3) = (a,w _1 /3) we obtain that r appears as the reflection 



in some root 7 on the root/weight side of the picture: rx = x — (7, x)6lq. The condition (21) 
is equivalent to 

(7, x) (A + p, do) S Z, for all x £ m^P ■ 

Note that (7, x) £ M + fe Z since (7, P) C Z, so we only need (A + p, do) £ (M + \i)7L. 
The simplest case is to look for (A + p, —do) = M + h, that is, 

n n 

J^iXj + l)( Uj , -do) = J2(\j + l)fhj = M + h, 
3=1 i=i 

or equivalently, 

n 

(22) ^ \jfhj = M + 1 . 

3=1 



All solutions A = (Ai, . . . , A n ) to (22), where the Xj £ Z-°, lead to ^-functions S\+ p that 



are zero at all EFOs of Ad-order M + h in the interior of the fundamental domain. 

6. Introducing the polynomial functions 

Following [T5], assign variables Xj to the characters for weights uij. Thus we have polyno- 
mial variables 

(23) Xi, X 2 ,...,X n , where Xj := Xui {*) , x £ F° . 
With these we can introduce the domain 

(24) := {(Xi(x), . . .,X n (x)) :x£F°}cC n . 

We shall soon see that this is actually an open subset of a real n-dimensional space and 
eventually it will be the natural domain of the functions real- valued functions of the variables 
Xi, . . . , X n that we wish to study. 

Define the Tn-degrea 5 ^ of the variables by assigning degree fhj to Xj. 



variable 


Xi 


x 2 




X n 


s P 


K 


m-degree 


rhi 


7T12 




rh n 


h-1 


2/i-2 



Table 1 . The ?7z-degrees of the polynomial variables X\ , . . . , X n and of the 
functions S p and K. The Coxeter number h = 1 + mi + • • • + rh n . 



This might be more properly called the rh-degree, but this seems a bit cumbersome. 
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Then the monomials X 1 1 . . . X^ n of m-degree < M are those satisfying 

(25) Aimi H h X n m n < M 

where Ai > 0, . . . , X n > 0. Although the marks and co-marks are not necessarily identical, 
Fig. [11 they are at worst simply permutations of each other. Thus (25) has the same number 



of solutions (Ai, . . . , A n ) as we saw before, namely \F^ I+h \ = \Fm\- The constant polynomials 
are those of m-degree 0. 

In keeping with this notation, we will say also that A = (Ai, . . . , X n ) has m-degree equal to 

(26) Aimi H h A„m n = (A, -do) • 

Theorem 6.1. 

The number of monomials C[X\, . . . , X n ] of m-degree < M is equal to the number of regular 
EFOs of Ad-order M + h in the fundamental chamber. Each of the regular EFOs of Ad-order 
M -\-h in the fundamental chamber is a common zero of all the S -functions S\+ P and all the 
character functions x\ f or which A has m-degree equal to M + 1. □ 

The trick that we have used above of using the internal reflective anti-symmetry to construct 
common zeros is taken from [8] . It is remarkable that in the case of type A n root lattices it 
actually finds all the common zeros. The proof of this makes essential use of the fact that 
the new variables X\, . . . , X n are all of degree 1, something that is true only for type A n . In 
fact our example of type G2 in ^9] indicates that this result does not hold in general. 

The 'smallest' S-function is the one defined by the strictly dominant weight of lowest m- 
degree, namely p of (171) with m-degree h — 1. Writing S p in its well-known form (p|, we note 
that ~S~ P = (-1)I A +IS P . Thus 

I c |2 _ c ~q _ /_-i\|A + | c2 

\Op\ — >->pDp — ^ L) Op , 

and we note that this function is positive on all of F° and vanishes on its boundary. S P S P is a 
VF-invariant function and so is expressible as a polynomial in the basic characters Xuij = Xj. 
We then have the corresponding strictly positive function K on £1: 

(27) K(X 1 ,...,X n ) = S p Sp(x), xeF° 
We note here that if n G P" 1 

S,, = Y^ (-l) l(w) e- Wfl =- ^2 (-l) l ( w) (-l) l{w °'»> ) e- ww °'*> fl = (-1)I A +I 5-t,.,.,.,./, • 

where w opp is the opposite involution in W (since w opp is a product of |A + | reflections). 

The opposite involution interchanges the positive roots (resp. positive coroots) with the 
negative ones, and, since fi is dominant so is —w opp fi. Of course w opp is not always simply 
the negation operator, see Tab. [2] Still, it does simply change the sign of the highest positive 
root (resp. highest coroot). Thus we have the important little equation 

(28) m-deg (/i) = {fx, -d ) = {w opp fi, -w opp a ) = {w opp fi, d ) = m-deg (-w opp jj) . 

This is useful because it means that (—ly A+ 'Sp and X/j. are J us t another 5- function and 
another group character respectively, and the highest weight involved in each case has the 
same m-degree as before conjugation. In particular conjugation of the characters \w can 
at worst permute some of them, say Xu>j *-* Xo-(w) by some permutation a of order 2 of the 
indices {1, . . . , n}. Table ^ shows what happens in the cases when a is not just the identity 
permutation. 



i) • • ■ iX n ) - 


: bpbp[x), : 


then 






E<- 


_!)»(« 


)(—iy( w op P ) e 


wew 
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If we let 

(29) m = {z = (z 1 ,...,z n ) :z = (z ct(1) , . . . , z^ n ))} , 

then 9\ is a real space of dimension n and 0, C IK is an ra-dimensional subdomain, as follows 



from the non- vanishing of the Jacobian S p on F° (31). This is the space on which we shall 
think of Xi, . . . , X n as real variables. 



X 


Xt. 


x 2 




X n -i 


x„ 




X 


X n 


X n -i 




x 2 


Xi 








X 


Xi 


x 2 




Xzn— 1 


X 2n 


X 2n +i 


X 


X\ X2 




X 2n -\ 


X 2n +i 


x 2n 



A„ (n > 1) 
D 2n+ i (n > 1) 

Table 2. Correspondence of the variables Xj to Xj, (J = 1, . . . , n) produced 
by the action of —w opp . In all other cases Xj = Xj. 



X 


x t 


X 2 Xs 


Xa. Xg 


X 6 


X 


x 5 


x 4 


Xz X 2 


Xi 


x a 



Remark 1. As we have seen, conjugation actually permutes some of the basic variables 
Xj. We shall use the overline symbol to indicate this form of conjugation. Thus one 
should understand the conjugation symbol as having this dual meaning of actual complex 
conjugation when the Xj are treated as functions on T and as the permutation a when 
treated as the coordinate variables of D\. Thus we shall write c Xj 1 . . . Xj r (where c £ 
C) to mean cXj 1 ...Xj r , understanding that the Xj has this dual meaning. For a poly- 

normal g(X^ . . , X n ) = £ c ju ..., jn X? ...Xfr, g(X 1 ,...,X n ) := £ cJ^^X? 1 . . . X~ n jn = 
g(X[ h ...JQ 3n ). 

Notice that since 



XjXj^K i > Xu>jXu> k bpbp — b w -J r pb UJk J r p 



and 



(ujj + p, -a ) = rhj + h— 1 and (-w opp (oj k + p), -a ) = m k + h - 1 , 
we understand that K has m-degree equal to 2h — 2. 



7. The integration formula 
We wish to study weighted integrals of the form 



f fgK l / 2 dX 1 ...dX n , 
Jn 

where /, g are functions of the variables X±, . . . , X n defined on Q,. These are related back to 



t (more specifically to F°) and the torus T via the defining equations (23). 
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7.1. The key integration formula. 

Natural variables for t are x = (xi, . . . ,x n ) = Yl 1 j=i x j^j where the Xj run over [0, 1) x 
•■■x[0,l). 

The derivation D^ on C[P], Da :j e x = (A,dy)e A is, when C[P] is treated as an algebra of 
functions on T, the mapping 

(30) = (A,a 7 )e< A ' 2 ™> = — — e^ 2mx K 

2tti dxj 



Using Prop. 3.3 we then see that the Jacobian of the transformation of the variables x to 
variables X is 

(31) \(2mTS p (x)\ = (27rT\S p (x)\. 

Thus from the definition of K we have 

f fgK l ' 2 dX l ...dX n = ( f(X 1} . . .,X n )g(X 1 ,...,X n )K 1 /\X 1 , ...,X n ) dX 1 ...dX n 
Jn Jn 



= ( 27r ) n / f(Xui (»),■••, Xu„{x)) g{Xui (»), • • • i Xco n (x)) S p (x) S p (x) dx x ... dx n 

J F° 

for all functions /, g are in the variables X±, . . . , X n on Q. 

Theorem 7.1. 

Let M be a positive integer. Then for all polynomials f,g€ C[Xi, . . . , X n ] with m-deg (/) < 
M + 1 and m-deg (g) < M we have 



[ fgK l / 2 dX l ...dX n 
Jn 



(32) 



( 27r ) n / /(Xwi(s), • • • , Xoj„ (x)) g(Xoj! (x),-.-, Xoj„{x)) S p (x) Sp(x) dxx... dx n 



1 / 2vr 



c G \M + h 



\ n 

) Yl /(^i( x )>--- 



XUn(B))9(Xui(B)> ■ ■ -,Xun(x)) S p (x)Sp(x) . 






This theorem is proved in £7.5 



7.2. The cubature formula. 

For any function / defined on Q,, let / be defined by f(x) = f(Xuiii x ), ■ ■ ■ ,Xui„( x ))- 

Corollary 7.2. 

Let M be a non-negative integer. Then for all polynomials f £ C[X\, ... ,X n ] with m-deg (/) < 
2M + 1 we have 



(33) 



/ /^«i-«- = 4(^) B E /»*(»)■ 



T cpo 
xt ' r Af+h 



Equation (33) is the cubature formula. The points x E F^f+h = ( M+h ) ^ -^° are 
the interpolation points. These points are common zeros of the character functions of the 
m-deqree M + 1. The coefficients of the interpolation are the values of K = \S P \ 2 at the 
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interpolation points. The Corollary is a direct consequence of Theorem 7.1 since every poly 



nomial of m-degree less than or equal to 2M + 1 can be written as a linear combination of 



polynomials of the form fg appearing in ( 32 ) . 
We first prove a key lemma. 

7.3. Separation lemma. 

Lemma 7.3. 

Ifcj) = J2] =1 4>jUj G P + and 4*^0, and ifm-deg{4>) < 2(M + h), then (f>(£(M + h)Q. 

Proof. Suppose by way of contradiction that <p E (M + h)Q. We have 

< (4>,-a ) < 2(M + h) 

from our assumption on the m-degree of eft and since ((f), —do) = Yl <A? ™j > 0- However, 
since (Q, — do) C Z, our assumption on cp forces J2 <pj rhj = ((/>, —do) € (M + h) Z and hence 

In the same way, applying </> to each simple coroot dj in turn, we obtain </>j = {cp, a.j) = 
(M + h) dj for some a^ G Z-°, j = 1, . . . , n. Thus M + h = J2 4>j rhj = (M + /i) X]?=i a j ^j- 
This implies that exactly one Oj 7^ and for this j, aj = 1 and rhj = 1. Thus = (M+/i)wj G 
(M + /i) Q, i.e. ojj G Q. 

In fact this can't happen. One way to see this is to use a well-known fact that the non- 
trivial elements of the centre of the simply-connected simple Lie group G are given by the 
elements ojj over the j > which correspond the places where rrij = 1 (these are certain 
vertices, different from the vertex 0, of the fundamental chamber. Of course these elements 
Uj $. Q for such an element would be a representative of the identity element. Our case here 
is the same situation except it is for the simply-connected simple Lie group G based on the 
dual root system: we have a j for which rhj = 1, and hence ujj ^ Q. 

Thus cf) = (M + h)ujj cannot lie in (M + h)Q. □ 

7.4. Weyl integral formula and its consequences. 

We recall here the Weyl integral formula [20J, 



1 1 -...,2 



^^0=7^7, / F\S P \ z d6 T 
G \W\ h 

for all class functions J- (functions which are invariant on conjugacy classes in G). Here the 
measures are normalized Haar measure on G and T respectively and the function T is simply 
being restricted to the maximal torus T in the second integral. In particular characters are 
class functions, and so for the irreducible characters X\->X P of the irreducible representations 
of G of highest weights A, [i respectively, we have from the standard orthogonality relations 

m- ' 

(34) 5 x ,n = / X\Xfide G = — / X\Xv\S P \ 2 d8 T . 

Here we are using the usual Kronecker delta. 



We do not need all of (34), only the equality of the left and right hand sides; and that fact 



is not hard to see. We have 

XxS pX r t tS p = S x+p S^ p = Y, ^l-iy^l-ll'We^^+^He-^W-*)) 
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by Weyl's character formul^ The integral over T of e ^M*+p)> a! >)- 2 * < («(H-p)»*>) ig unlesg 
w(X + p) — v(p + p) = 0, in which case it integrates to 1. Since A + p and p + p are strictly 
dominant, this happens only if A = p and w = v. If indeed X = p then there are exactly \W\ 
times when w(X + p) — v(p + p) = 0, and we see that the right hand side of (34) is 5\ ifl . 



We also wish to recall a result from discrete Fourier analysis [13] • There it is proved that, 
using the notation established above, 

(35) S x+p S IM+p d6 T = — — V S x+P {x)S p+p (x) = 5 x ,p , 

Jf c G (M+h) n J-f 

xt - r M+h 

as long as when A / p, the points of F^ +h , can separate all the weights appearing in the 
W-orbit of A + p from all those appearing in the PF-orbit of p + p. Explicitly, separation 
means that it never happens that w(X + p) — v(p + p) takes integer values on all the points of 
M \ h P, or equivalently, it never happens that w(X + p) — v(p + p) 6 (M + h) Q, except when 
X = p. 



This proof of the last equality in (35) is actually a straightforward thing to see. First note 
that Sx+pS^+p is IF- invariant, and hence its integral over all of T is \W\ times its integral 
over F. There is no need to worry about the boundary of F which has measure and in any 
case the function takes the value on all of its boundary. In the same way, the sum over F^ +h 
can be extended by the operation of the Weyl group to obtain a full set of representatives 
of M \ h P/Q, noting again that since the function is on the boundaries of the chambers, 
adding in the extra boundary elements that may appear in M \_ h P/Q makes no difference to 
the sum. This again increases the value of the sum by \W\ but has the benefit of turning the 
sum into a sum over a group. Then usual considerations of sums of exponentials over a group 
give the desired orthogonality as long as the separation condition is satisfied. The order of 
jrrv P '/Q is cg (M + h) n , which explains the factor outside the sum part of the formula. 

7.5. Demonstration of Theorem 17.11 



With these facts, we now prove the integral formula (32): 



Proof. Because of linearity, it is enough to prove ( 32 ) when /, g are monomials of the form 

X^\ Xul ■ ■ ■ Xul > where ^i m i + ^2"i2 H \- v n m n < N , 

where N = M + 1 for / and N = M for g. 

Since XwiXoS • • • X^™ decomposes into a linear combination of characters ^ a AXAiw 1 H hA n w n , 

where Xj < Uj for all j, and a u = 1 (i.e. the exponent v = v\uj\ + z^W2 + • • • + v n LLi n appears 
in the sum with multiplicity 1), we need only prove the Theorem for f(Xi, . . . ,X n ) = Xx = 
XAiwi+-+A»w„ and g(Xi, . . .,X n ) = Xp = Xtnui+-+HnUn- Notice that it is quite possible for 
X0) which is the constant function It, to appear here. 



We have from (34) and (35) 



(36) 8x,fjt = TTT^j- / X\Xn\Sp\ dOf = —=- / Sx+pS^+p d9j = / Sx+pSp+pdQj 

1 



cg(M + h)' 



/ j >->A+p\^/^n+p 



b\+p{X)bp+p{X) , 



In [20] Weyl's character formula is derived from the integral formula, but algebraists usually use an algebraic 
proof of the character formula. 
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(which is exactly what we have to prove) as long as when A 7^ p there are no pairs w, v £ W 
for which w(X + p){x) — v(p + p)(x) G Z for all x G M \_ h P, or, as we pointed out above, 

w(X + p)-v(p + p)e(M + h)Q. 

We will now show this cannot happen. 

Consider the weights w(X + p), w E W. These weights are all of the form A + p — /3 where 
/? is a sum of positive roots (including the case when it is the empty sum, 0). Now —a is the 
highest co-root and so is actually a dominant co-weight. This gives us that (/3, — a) > and 
so 

m-deg (w(X + p)) < m-deg (A + p) < M + h . 

The lowest weight in the VF-orbit of A + p is w opp (X + p) and its m-degree is the negative 
of the m-degree of A + p as we saw above. All the other weights in its orbit are of the form 
Wopp(X + p) + (3 for some sum (3 of positive roots and this then gives us 

— (M + h) < —m-deg (A + p) < m-deg (w(X + p)). 

In short 

-(M + h) < -m-deg (A + p) < m-deg (w(X + p)) < m-deg (A + p) < M + h . 

Exactly the same holds for v(p + p) except that the inequalities are now strict since the degree 
of g is at most M. Combining, we obtain 

-2(M + h) < w(X + p) - v(n + p) < 2(M + h) 

for all w,v £ W. 

Now for any choice of w, v, each element in the W-orbit of w(X + p) — v(p + p) is another 
element of the same form, and so its m-degree is constrained in the same way. Thus if there 
is a pair w, v for which w(X + p) — v (p + p) G (M + h)Q then, since (M + h)Q is VT-invariant, 
we can assume that 4> '■= w(X + p) — v(p + p) E P + n (M + h) Q, i.e. it is dominant. But then 
if 7^ it contradicts Lemma |7.3| It follows that (j) = 0, and then due to the fact that X + p 
and p + p are strictly dominant, we have X = p. 

This proves the separation condition holds, and finishes the proof of the Theorem. □ 

7.6. Duality. 

The essence of the cubature formula of Corollary 7.2 is the duality between dominant 
weights of m-degree not exceeding M and the elements of finite order M + h arising from the 
fundamental region F. If we use the fact that the m-degree of p is h — 1 and the fact that 
any regular EFO in F can be expressed as p + p (or x^+p in T), where p is co-dominant, then 
the cubature matrix is : 

(M A A ) = (Sx+pixji+p)) , 
where A, p run over all solutions to the equations 

^ThjiXj + l) = M + h 
^mjipj + 1) = M + h. 

Since the marks and co-marks are simply permutations of each other, the symmetry of this 
pairing is completely manifest: if we order the co-ordinates on each side to take account of 
this permutation, the solutions to the two equations look identical and the matrix becomes 
symmetric. Moreover, formally the situation is the same for G and the corresponding group 
G, with the roles of characters and elements of finite order interchanged. 
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8. Approximating functions on 

To simplify notation and the following discussion we introduce an inner product on the 
space L^(fi) of all complex- valued functions / on $1 for which J n \f\ 2 K 1 ' 2 < oo. It is defined 
by 

(37) 

(f,9)i<-= f fgK 1 ' 2 = f f(X x , . . .,X n )g{X u ...,X n ) K 1/2 (X l5 ...,X n ) dX 1 ...dX n . 
Jn Jn 

Note here that the definition of conjugation is made in terms of Remark [TJ 

Since K 1 ' 2 is continuous and strictly positive on O, (/, /)k > with equality if and only 
if / is zero almost everywhere in fi (relative to Lebesgue measure). Relative to this inner 
product L 2 K (Q) is a Hilbert space. We call (/, f)ji the L^-norm of /. 

The results of (36) show that the polynomials defined by 

X x := X x(x), xeF°, 
where A = ^ AjWj runs through the set of dominant weights P + , form an orthonormal set: 

(X\,Xfj)K = o~\,n ■ 

These polynomials are just those that result from expanding xx as a polynomial (of degree 
equal to the m-degree |A| m of A) in the fundamental characters Xj = Xuj- 

The set {X\} of actually forms a Hilbert basis for L 2 K {Si) (the main point being that 
they actually span the entire space). This can be seen by relating functions / on Q back to 
functions /on F° through f(x) := f(X 1 ,...,X n ) = f{Xu>i(x), . . . ,X w „(x)). Then 

(38) oo> f \f\ 2 K 1 / 2 = f \f\ 2 (x)SJx)SJxjdx. 



n Jf° 



Formally / exists only on F° , but we can extend it to a function on all of T by PF-symmetry 
if necessary. This might make it a bit easier to relate to £J3j which we are now going to use. 
We have already seen the integral on the right hand side of (38) in (14), and according to 



(12) and (13) we have 

/ = J2 bx Xx = ^2 &A Xx ' 

A A 

with 

bx = TWT [ fXx\Sp\ 2 d6 T = [ fxx\S P \ 2 d6 T = [ fX~ x K x l 2 = (f,X x ) K , 
I w I Jf Jf° Jn 

(see 0). 

Thus for each / £ L 2 K (&) we have its expansion 

f^^2(f,X x ) K X x , 

A 

where the — means equality in the sense of equality Lebesgue-almost-everywhere. 
The truncated sums 

/ ] (f, X\)kX x , 

\\\m<M 

where |A| m stands for the m-degree of A, are polynomials of m-degree at most M in the 
variables X±,..., X n . 
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8.1. Optimality of the approximation. 

We now show that these polynomials are the best possible approximations for / in terms 
of the L^-norm. This result is in fact a natural consequence that one would expect from the 
situation that we have created here, since it essentially relates back to the Fourier analysis of 
T. See also EH. 



Proposition 8.1. 

Let f G L 2 K {Q). Amongst all polynomials p(X±, . . . ,X n ) of m- degree less than or equal 
to M, the polynomial X^iai <m (/>-^a)k X\ is the best approximation to f relative to the 



l? K -norm. 



Proof. Let p = X/A| m <M^A-^A be any polynomial. Let a\ := (f,X\)x for all A and set 

Q. '■= E|A| m <A/«A^A- 

(/ ~P,f ~P)k = (f,f)K ~ (f,p)K - (pJ)k + (p,p)k 

= (/,/>k-J>a&a~-J>aoa + 5>a| 2 



AAA 



(/^k-EH'+EM'-E^-E^+EN' 

A A A A A 

(/ - q, f - q)K + E I 6a ~ flA ' 2 - (f ~ 9 ' f ~ qS > K 



A 

with equality if and only if 6a = a a for all A. □ 

9. Example: A cubature formula for G 2 

We consider an example where the Lie group is the exceptional simple Lie group of type 
G2 and the value of M is chosen to be M = 8. Before we really start, let us recall pertinent 
information about G2 and related objects. Then we can start to solve the problem of our 
example. There are three steps to it. First the S- functions of m-degree M + 1 = 9 are to be 
found (or more specifically the highest weights that label them) . Then their common zeros of 
Ad-order M + h = 14 in the interior of F are to be determined, and then only the cubature 
formula (matrix) can be written down. 

9.1. Pertinent data about Gi- 

From the G2 diagram in Figure [T] we find the following: 

(ai,ai) : (a 2 ,a 2 ) 
(&i,&i) : (02,02) 
— ao = 2«i + 3«2 5 
marks 2, 3, 
Coxeter number h = 6 . 
Therefore we have 

ct\ = 2oj\ - 3w 2 , a 2 = -ojx + 2ui 2 , 0J\ = 2a± + 3a 2 , w 2 = «i + 2a 2 , 

&i = 2a)i — w 2 , a 2 = — 3tDi + 2cD 2 , Co\ = 2&i + a 2 , w 2 = 2>a\ + 2a 2 . 

The link between the bases is the Z-duality requirement 

(aj,LJk) = (a.j,u k ) = 5 jk , j, k = 1, 2 . 



3:1, 


A=( 2 r 2 3 ); 


1:3, 


A = (-3 2 ) 1 




—do = 3di + 2a 2 




co-marks 3, 2 
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Figure 2. A schematic view of the co-root system of G2- The shaded triangle 
is the fundamental region F. The dotted lines are the mirrors which define 
its boundaries, the reflections in which generate the affine Weyl group. The 
action of the affine Weyl group on F tiles the plane. A few tiles of this tiling 
are shown. Filled (resp, open) squares are the short (resp. long) co-roots of 
G 2 . 



The set of positive roots A + and their half sum p are: 

A + = {ai, 02, a.\ + «2, Oi\ + 2a 2 , «i + 3a2, 2«i + 3a2J, p = 3ai + 5a 2 = loi + uj 2 ■ 
The fundamental domain of G 2 is the convex hull of its three vertices {0, g&i, 5^2}- 

9.2. Finding S'-functions of m-degree 9. 

The ?TT,-degree of a G 2 weight A = {Ai, A2} = Aiwi + A2W2 is calculated as miAi + 7712X2 = 
3Ai + 2A2- We find first all the weights with m-degree < 9. They are the following, 

{0,0}, {0,1}, {0,2}, {0,3}, {0,4}, {1,0}, {1,1}, {1,2}, {1,3}, {2,0}, {2,1}, {3,0}. 

Among these there are just two weights with m-degree equal to 9, namely {1,3} and {3, 0}. 
Thus there are only two S- functions we need to consider, Si^(x) and Sso(x). 



9.3. Finding common zeros of the 5-functions of m-degree 9. 

The interpolation points are the points of the set F°± of EFO's of Ad-order M + h = 14 
which are in the interior of F. They must satisfy the sum rule 14 = sq + 2si + 3s2- We write 
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them in the in Kac coordinates [sq, s±, 8%] as well as in the w-basis: 

[9,1,1] = (ti,t3), [7,2,l] = (i,i), [5,3,1] = &,u), [3,4, 1] = (f , i), 
(39) [1,5,1] = (£,£), [6,1,2] = (i,l), [2, 1,1] = (1,1), [2,3,2] = (^,1), 

[3,l,3] = (i,^), [1,2,3] = (£,&) 
The strict adjoint order of all the EFO's but one is 14. It is 7 for l&i + ^U2- Other EFO's of 



adjoint order 14 are not shown in (39) because they are on the boundary of F. The boundary 
points are easily discarded because at least one of their Kac coordinates has to be zero. At 
these points the 5-functions Si^ and S^o simultaneously vanish. 

9.4. Cubature formula. 

Underlying the cubature formula (for M = 8 in G2) is the fact that the matrix 

X = X® = (Sx+pix/,)) 

where the A run over the 10 weights of m-degree at most 8 and the fi run over the 10 regular 
EFOs of Ad-or&ev 14, satisfies j^—^oXX T = iioxio- Here is the matrix X^> with entires 
rounded to 4 figures. 

-0.604 -2.714 -5.494 -6.098 -2.714 -3.384 -7. -6.098 -3.384 -2.11 

-2.714 -7.604 -6.098 1.506 2.714 -4.89 0. 6.098 3.384 3.384 

-5.494 -6.098 2.11 -3.384 -6.098 2.714 7. -3.384 2.714 -0.604 

-6.098 1.506 -3.384 -4.89 6.098 7.604 0. 3.384 -2.714 -2.714 

-2.714 2.714 -6.098 6.098 -7.604 3.384 0. 1.506 -4.89 3.384 

-3.384 -4.89 2.714 7.604 3.384 -1.506 0. -2.714 -6.098 -6.098 

-7. 0. 7. 0. 0. 0. -7. 0. 0. 7. 

-6.098 6.098 -3.384 3.384 1.506 -2.714 0. -4.89 7.604 -2.714 

-3.384 3.384 2.714 -2.714 -4.89 -6.098 0. 7.604 -1.506 -6.098 

-2.11 3.384 -0.604 -2.714 3.384 -6.098 7. -2.714 -6.098 5.494 

Direct computation shows that X^(X^) T = 14 2 /i 0x io- 

The values of the weighting K function at the interpolation points are given by 

0.364666, 7.36467, 30.1836, 37.1836, 7.36467, 11.4517, 49., 37.1836, 11.4517, 4.45175 

The image of the fundamental domain F under the mapping x \— > (Xi(x),X2(x)) in IR 2 
along with the images of the 8 EFOs that we have used in this example is shown in Fig. |3J 
The way in which EFOs fill out the region il is made clearer by the Fig. [4] which shows the 
distribution of the EFOs of Ad-order 106: 

Examination of the graphs of the two functions £1,3 and Ss,o reveals that they have (at 
least) 2 other common zeros in F. These do not seem to be related to EFOs. 
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Figure 3. The region f2 along with the 10 regular EFOs of the example. 




Figure 4. The region Q along with the 884 regular EFOs of Ad-order 106. 
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